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ABSTRACT 


The cell contact with a solid substrate is of great importance in the study 
of developmental biology. The interaction between a cell-membrane and a solid 
support can be modeled by considering the stability of the fluid film that is held 
in between them. This fluid film is a viscous layer which is subjected to forces 
like van der Waals interactions and the polar (acid-base) interactions. In the 
present study, two cases are considered for the stability analysis of the film, - i) 
effect of bending resistance offered by the cell membrane during its deformation, 
and ii) effect of aggregation of receptor molecules on the membrane surface. The 
movement of the membrane-film interface is described using a hydrodynamic ap- 
proach. A systematic long-wave analysis is done to derive the nonlinear evolution 
equation governing the position of the said interface. For the second case, this is 
coupled with the evolution of the receptor concentration (based on the principle 
of conservation of mass). The non-linear evolution equations have been solved 
using appropriate numerical techniques. The results with different parameters 
are presented for the two cases. 
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NOMENCLATURE 
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E nondimensional amplitude of perturbation 

f friction factor 

G free energy 

h dimensional film thickness 

H nondimensional film thickness 

k dimensional wave number of perturbation 

K nondimensional wave number of perturbation 

k c dimensional critical wave number of perturbation 

K c nondimensional critical wave number of perturbation 

k m dimensional dominant wave number of perturbation 

K m nondimensional dominant wave number of perturbation 

k mi dimensional dominant wave number considering 
interfacial tension only 


5 



k m b dimensional dominant wave number considering 
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l nondimensional correlation length for the polar forces 

M moment produced by bending 

p film pressure 

P nondimensional ratio of polar to apolar spreading coefficients 
i?i, i ? 2 radii of curvature 

R\ m ratio of the dominent wavelengths for the case with 
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S LW spreading coefficient for the apolar force 

S AB spreading coefficient for the polar force 

t dimensional time 

T nondimensional time 

Ti nondimensional time of rupture obtained from linear stability analysis 
x dimensional space coordinate 

X nondimensional space coordinate 
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e dimensional amplitude of perturbation 
7 interfacial tension 
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A dimensional wavelength of perturbation 
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interfacial tension only 

A nondimensional wavelength of perturbation 
H viscosity of the film fluid 

v kinematic viscosity of the film fluid 
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Chapter 1 


INTRODUCTION AND 
LITERATURE SURVEY 


In the progress of developmental biology, the determination of the developmen- 
tal fate of a cell as a result of interactions with its environment has widespread 
applications. The environmental interaction is often caused by the specific direct 
influence of a solid support ( extra-cellular matrix ) on the cell. This type of in- 
teraction is the basic to such prosseses as histogenesis, wound healing, metastatis, 
neuronal organisation etc (1). Cell contacts made as a result of these processes 
may play important role in cell physiology. They are involved in the cell locomo- 
tion as well as in the extensive spreading and immobilization of the cell. They 
also contribute to the acquisition of the cell polarity, transduction of the tran- 
scellular mechanical forces, formation of the permeability barriers, control of cell 
growth etc (2). All of these point to the necessity of the study of the basic physics 
for adhesion of cells. 

Several workers in the field of bio-physics have visualised this adhesion phe- 
nomenon in the light of thin film ( the liquid films being thin enough so that 
their thermodynamics and trasport behaviours can depend on their thickness 
which generally ranges from 1 to 100 nm) concept of interfacial science. These 
two fields are combined to model the cell-support adhesion by considering the 



stability of the liquid film in between the cell membrane and the support. The 
stability of the thin film has been analysed by employing linear theory considering 
only the van der Waals force as the body force (3). The nonlinear approach using 
long wave analysis was subsequently adopted (4). This is then further modified 
through incorporation of the polar forces and other specific interactions. This 
concept has been applied to study the cell-support interactions. The prediction 
of the most probable wave length for such adhesion (5), the tear film rupture anal- 
ysis (6), the study of presence of receptors on the cell membrane using Marangoni 
effect in a non-linear approach (7) etc. can be cited as examples. 

In the present work two specific characteristics of some cell-membranes are 
taken into consideration for the stability analysis of the film, namely the effect of 
resistance offered by the membrane due to the bending moments developed in it 
during deformation and due to interfacial tension (interfacial interactions), and 
the role of receptors present on the membrane surface. Bending moments are 
generally considered to stabilize the process of cell adhesion, providing a stabiliz- 
ing effect. This effect is studied together with the interfacial tension, for different 
parameters of the system ( like the body force term, mean film thickness etc.) 
to predict the actual cell-support contact or the deformed equlibrium structure 
in absence of cell-support contact. There are many well documented systems in 
which the presence of some specific molecules ( lectins, SAMs, antibodies etc) is 
required for the cell-support adhesion to occur (8). It has been proposed that 
these molecules can overcome the repulsive barrier for which cell adhesion is not 
possible in case of bare membrane (9). It has been considered that they can form 
specific bonds at a short distance near the position of actual contact. In this 
work this phenomenon is studied in the light of effective free energy due to the 
presence of these receptors, for various system parameters. 
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Chapter 2 


EFFECT OF INTERFACIAL 
INTERACTIONS IN 
CELL-SUBSTRATE 
ADHESION 

2.1 INTRODUCTION 

In the process of cell adhesion, the stabilising effect provided by the interfacial 
tension and the bending resistance is important. In order to account for the elastic 
rigidity of the membrane, the effect of bending moment has to be considered. 
Bending moment resists the formation of curvature on the surface (10). The effect 
is especially pronounced for the cases involving significant curvature changes. The 
present work involves the study of this influence of bending moment during the 
adhesion of a cell membrane on a solid support along with interfacial tension. It 
considers different types of cell membrane, support and the fluid film sandwitched 
in between them. This leads to the study of the system with different cases of the 
interaction potential considered to describe the body force, and varying interfacial 
tension. 
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2.2 CONCEPT OF BENDING MOMENT IN 
CELL MEMBRANE 

Let us consider the membranes with connected bilayer (by van der Waals contact). 
In this case the bending moment at any point on the cell membrane and film 
interface can be related to the local curvature (10). If the connected bilayer 
is initially flat and then bent, a bending moment is produced by the effect of 
expansion of the upper layer and compression in the lower layer causing a small 
area change. This moment can be expressed as (10): 

< 2 - ! > 

where B is the bending rigidity constant, and is a function of the membrane 
properties like compressibilities of the layers and the distance between them. Ri 
and R 2 are the local radii of curvature. The force due to this bending moment 
can be obtained as V 2 M. Now this force is incorporated in order to balance the 
normal stress equilibrium condition at the cell-film interface given as follows: 



Here A p is the pressure difference at any point on the membrane-film inter- 
face, and 7 is the interfacial tension of that interface. 

2.3 FORMULATION 

A cell separated by a thin liquid film from a solid support is considered. These 
are separated by an interaction region which can be visualised as a, non-thinning 
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film of thickness h , bounded laterally in the contact region with co-ordinate 
x = 0 and x = x. The two interfaces are at 2 = 0 (support and film interface)and 
2 — h, (cell-membrane and film interface) (fig. 2.1). 

Starting from a reference state with a plane film of mean thickness ho the 
movement of the fluid inside the film is described by the Navier-Stokes equa- 
tion for anisotropic, incompressible, non- draining and viscous fluid of viscosity p 
and density p. The force balance is between apolar (attractive) and polar forces 
(attractive or repulsive). Dynamic boundary conditions are considered at the 
cell-film interface while the velocities are zero at the support-film interface.The 
boundary condition on the normal stress involves the product of the surface ten- 
sion and the curvature effect, and bending moment of the cell membarane. The 
shear stress on the cell-film interface vanishes due to constant surface tension. 

A long wave analysis (4) for this thin film (with the condition |<1) leads 

au 

to the evolution equation for the cell membrgfie (2 = h) as follows: 

h ‘-W^ h TJ» + » ] = 0 (2 ' 3) 


where the p and 4> are given by (using small slope approximation): 




p — y hxx T ^ h 
2 S LW cP 0 S AB 


h* 


lo 


.do — h. 

ex P(— 7 ) 

‘0 


(2.4) 

(2.5) 
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Here 7 is interfacial tension at cell-film interface , B is bending moment of the 
cell membrane, S is spreading coeffcient of the film on the support in presence 
of the cell, LW refers to apolar and AB to polar interactions, d 0 (— 0.158 nm) 
(11,12) is the equilibrium distance for the film where the born repulsion may be 
replaced by a vertical rise in the potential to infinity (hard sphere approximation) 
and lo (— 0.6 nm) (11,12) is the correlation length for the polar forces. 


NON-DIMENSIONALISATION: 

The above equations can be non- dimen sionalised using the following non- 
dimensional parameters to retain mininmum number of parameters and for com- 
pact representation: 


nr ^ j ^0 , lo 

H = — ,d= — = — , 
ft 0 fto fto 


p = 


S AB 


S LW (6d*Py 

and non-dimensional space and time scalings, 


Y _ «£ r - — 
ho' hV 


where 


( 2 . 6 ) 


a = 


S LW <Ph 2 0 


B 


d 2 h 0 S LW a 2 

2v 2 p 


v = — 


(2.7) 


Now the evolution equation can be written in a form that involves least num- 
ber of parameters, namely, P, d , l and A which solely depend on surface properties 
(7, B , S LW , S AB ) and the mean film thickness h 0 : 
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H t = 


H 3 (-AHxxx + Hxxxxx) + sgn (S lw ) (iT 1 - P exp H 3 ) H x 

( 2 . 8 ) 

where the function s gn(S L ™' ) refers to the sign of S LW i.e, s gn(S LW ) = +1 if 
S LW > 0, and for S LW < 0 it is —1 and A is given by: 

A = ^ (2.9) 

(6 \S LW \Bd?)> 

The second and third term of the equation [2.8] are related to the stabilizing 
influence of the film interfacial tension and bending moment respectively. The 
remaining non-linear terms which account for the excess free energy (</>) of the 
film, may be stabilizing or destabilizing, depending on the signs of S LW and S AB . 
When S LW is negative Lifshitz-van der Waals forces are destabilizing and vice- 
versa. The same applies for the polar interactions based on the sign of S AB . The 
dynamics of the film is therefore dictated by the balance of nonlinear effect of 
the film tension and bending moment and the apolar and polar intermoleculer 
forces. The role of the non-dimensional parameter A is in studying the effect of 
the interfacial tension together with the effect of bending moment. 

2.4 LINEAR STABILITY ANALYSIS 

The initial (short time) evolution of the film profile is adequately described by 
linearized equation whenever the amplitude of the initial disturbance is much 
smaller than the mean film thickness. In order to get an initial estimate of the 
behaviour and the associate time scale of the film for infinitesimal perturbations, 


15 



linearization around the mean film thickness can be done. For the linear evolution 
equation [2.3] this mean film thickness can be taken at h = ho and for the 
nondimensional [2.8] one this will be H = 1, are done. Linearization leads to the 
following solutions respectively: 


h = ho + t sin kx exp (c ot) 

(2.10) 

H = 1 + EsmKXexp (f IT) 

(2.11) 


where e is the dimensional amplitude (e <C h 0 ) and E is defined as the ratio of 
e to ho- The dimensional growth coefficient to is related to the dimensional wave 
number, k as (4): 


u = 


3 fj, 


- h 3 0 ( 7 i 4 + Bk e ) - hlk^ (fa)," 


( 2 . 12 ) 


where {<j> h )h 0 refers to the derivative of 6 with respect to h, evaluated at h 0 . 
Similarly the non-dimensional growth coefficient ft can be expressed by the 
following equation: 

n = -AK 4 - K 6 - sgn (5 LW ) 


1 — P exp 


d- 1 


K‘ 


(2.13) 


where K is the nondimensional wave number. 

Now the critical wave number k c (0 < k < k c ) is defined for to = 0, since the 
surface instability can only grow if u> > 0 : 
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(2.14) 


k c 



-7 + \ff - 4 B 
25 


According to the linear theory, the fastest growing initial disturbance has a 
wave number k m (known as dominant wave number), obtained from the solution 
of |x = 0 , which gives: 


km. — 


N 


-7 + yf? - 3 5 

SB 


(2.15) 


Two asymptotic cases can be derived from the above equation, as shown 
below: 

For B -> 0 

= ( 2 . 16 ) 


\ 2 7 


and 7 — > 0 


k — 

“'rri'y — 




ho 


SB 


( 2 . 17 ) 


Similarly the expressions for the nondimensional critical and dominant wave 
numbers can be found from the following equations respectively: 



j A 2 — 4 sgn ( S LW ) 1 

(l — P exp | 

(¥)) 

N 2 


( 2 . 18 ) 
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Thus the non-dimensional wave length (in X coordinate) for the dominant 
mode of disturbance is given by A m = and its rate of growth (in T coordinate) 
is given by [2.13] Replacing K by K m . The dimensional case can be treated in a 
similar way. 

From linear theory one can predict a. hypothetical time of rupture (nondimen- 
sional) by setting H = 0 at sin (KX) = —1. This gives a linear estimate for the 
time of film rupture Tl, for a perturbation of wave number K , given by: 

Tl = 



The dimensional case can be treated in the similar way for evaluating domi- 
nant wave length and the linear time of rupture. 

2.5 NONLINEAR ANALYSIS 

Since the linear theory underestimates the nonlinear terms in the evolution equa- 
tion at the cell-film interface, it fails to predict either qualitatively or quantitavely 
the ultimate fate of the system. This necessitates the nonlinear analysis of the 
system. Here to understand the effect of the bending moment of the cell mem- 
brane on the evolution equation, the entire study is done in two parts, namely, 

1. The cases with no polar interactions (P = 0) 

2. The cases with polar interactions.(P ^ 0) 
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For all of the above cases the effect of interfacial tension is observed by varying 
the values of 7 at constant h 0: and the effect of mean film thickness is studied by 
the reverse way. The value of A is also varied as a lumped parameter for some 
cases. 

NUMERICAL METHOD: 

The sixth order nonlinear equation [2.8] may be solved with periodic boundary 
conditions over a wavelength, A = ^?,viz., 


'd'H' 


f d i H\ 


,(* = 0,1,...,5);0 < X < A 


( 2 . 21 ) 


\ax‘) x=0 \dx<) x ^ 

and a periodic initial condition can be used on the basis of initial perturbation 


as: 


H{ 0,X) = l+EcosI<X,(\E\ < 1) 


( 2 . 22 ) 


Cosine perturbation is taken as initial condition due to its symmetric nature 
about the mid-plane of the wavelength. This even reduces the computational 
time, since in this case the actual simulation for half the wavelength is sufficient, 
the other half being the mirror image of that. Fourier collocation technique (13- 
15) is used for discretization. The number of grid points required for accuracy 
is around 50 in the half wavelength. It is the most straight forward method for 
the problem with spatial periodicity. An automatic adjustment of the time step 
was employed by using the Gear algorithm for time integration. This becomes 
increasingly important near the time of rupture (16). The periodicity interval is 
first normalized to (0,2tt) by rescaling the X coordinate as: 

t) = KX, 0 < 7 < 2tt (2.23) 
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The discretization of the evolution equation and initial and boundary con- 
ditions are done in the following way (17). Eq. [2.8] , rewritten in 77 coordinate, 
is then exactly satisfied at N (even integer) number of Fourier collocation grid 
points (the end point is omitted due to periodicity). 


V j = 2tt (j - 1) /TV, j = 1,2, , N. (2.24) 


The collocation points are related to the discrete Fourier coefficients of the 
film thickness, a* 


N 

a k = (1/AO E H (Vi) exp"^, - (N + 1) /2 < K < {N + 1) /2 (2.25) 

i=i 

In view of this, Fourier collocation differentiation can be represented by a 
N x N skew-symmetric matrix, T>n with components 

(D w )„ = (1/2) (-1)'« cot [(>' -j) n/N] ,i/j (2 26) 

The higher order derivative matrices, D^, D^- and D^, etc. are simply 
obtained by n multiplications of the matrix D v , where n equals the order of 
derivative. 

The requirement that eq. [2.8] is exactly satisfied at each of the collocation 
points, gives a set of N equations. 


dH 

dT 


+ AK* [D„ {H 3 (DJtH)}] + K e [Djy {H 3 (DJ,H)}) 


[t„Hj (D jV H)] 
where H is a column vector of film thickness, 


(2.27) 
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H H (rj2, T ) , ,H( Vn ,T)). (2.28) 

H 3 is a N x N diagonal matrix with jth diagonal element equal to Hf ( r\j , T), 
and $h is iV x iV diagonal matrix with kth diagonal element equal to 

$hk = - Pexp {(d - H k ) //}] , (2.29) 

The N initial conditions are written as: 

H t -(?7 t -,0) = 1 + E cos (rji) , i = 1,2, ,iV (2.30) 

The time integration is done using the package D02EJF from NAG library 
and it is continued until the film goes to rupture (when at some localized point 
the minimum film thickness is falling below the cut-off thickness, d 0 ) or assumes 
some time independent profile. 

2.6 RESULTS AND DISCUSSION 

Based on the parameter values (11,12,18) of tables (2.1) and (2.2), the spreading 
coefficients S LW and S AB are calculated for different cases from the following 
equations, using the concept given in the appendix of (17): 

= 2 [V7^ - VtF] [fiF - \Ffi\ < 2 - 31 ) 

+ v^-i/7/j + /i7(v^ + V7+-v^) l (232) 

- ~ \Ac + 77 J 
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Table 2.1: Different components of 7 in mJ/m 2 for different cells and water 



'^rr 

7 + 

7" 

Erethrocytes 

35.2 

0.01 

46.2 

Lymphocytes 

27.7 

0 

32.2 

Corneal epithelial 

: 

27.4 

3.2 

56.3 

water 

21.8 

25.5 

25.5 


Table 2.2: Different components of 7 in mjjm 2 for different substrates 



, yLW 

7 + 

7" 

Agarose 

41 

0 

34 

Polystyrene 

42 

0 

1 

Cellulose acetate 

38 

0 

33 

teflon 

17 

0 

0 





where the suffixes c, /, and s denote cell, film, and support respectively. LW 
and AB refer to the apolar and polar components respectively. The typical value 
of the bending modulus B can be taken as 10~ 13 ergs (7) for all the cases described 
below. 

The relative importance of bending moment and interfacial tension will be 
clear if the dominant wavelengths for these two cases are compared. The value 
of dominant wavelength [\ m B — F%)j f° r the case with bending moment and no 
interfacial tension, is obtained by using eq. [2.17]. 


(2 - 33) 

Similarly, for the case with interfacial tension and no bending moment the 
dominant wavelength can be written as [2.16]: 

a ” 7 A(A) (2 ' 34) 

The ratio ( R \ m ) for these two dominant wavelengths, corresponding to the 
above two cases is given as: 


D A mB _ /3 B(-<j>h) hQ \* 

Rx ~ = - [ " 4 ?—) 

This can be expressed in terms of nondimensional quantities: 


(2.35) 


A 

' m y/A 


where: 


(2.36) 


Po = 


-!*»(»") (i- p «= x p(V)) 


(2.37) 
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and A is given by eq. [2.9]. 

For different sets of values of 7 and h 0 (at constant B ), the values of A m7 
are listed in table (2.3). Typical values of S LW and S AB are taken to be the 
same as those for a system with water film separated by polar cell-membrane 
and apolar solid support. It is clear from the table that for low values of 7, the 
effect of bending moment has to be included in the stability analysis of the film. 
Otherwise the length scale of instability would be much lower than the typical 
dimensions of the cells (0.01 to 1 /zm). Also the theory of long wave analysis is 
not found to be valid for these cases, since the wavelengths are smaller than the 
mean film thickness. However, ho remaining constant, the values of \ m B are of 
the same order as the cell dimensions (shown in table 2.4). Long wave analysis 
is also seen to be applicable in this case. On the other hand, for higher values of 
7 with high ho, the effect of bending moment may be neglected. The variation 
of R\ m with typical values of A and Pq is shown in fig. (2.2). The low values of 
A and high values of Po , characterised by low 7 and h 0 , correspond to the high 
values of R\ m . 

Now we can broadly classify the study into two parts. This is based on the 
type(s) of potential selected to describe the free energy of the system. 

2.6.1 CASE : 1 P = 0. 

In this case only apolar interactions are taken into account. This is the case when 
the film fluid and either the cell membrane or the support are apolar, i.e. both 
the polar components ( 7 + and 7“ ) for them are zero (e.g. hydrocarbon films in 
contact with lymphocyte cell). S LW has to be negative in this case to describe 
the attractive behaviour of the van der Waals force. 
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Table 2.3: Values of X mi for different 7 at different h 0 , S LW = 2 mJ/m 2 and 
S AB = —33mJ/m 2 . 


7 

(mJ/m 2 ) 

ho 

(nm) 

Am 7 

(H 

0.001 

2 

0.00013 

0.001 

8 

0.0256 

0.1 

2 

0.0013 

0.1 

8 

0.2 

1 

2 

0.0043 

1 

8 

0.8113 


Table 2.4: Values of A m g for different for different ho, S LW = 2 mj/m 2 and 
.S'" 4 ® = —33 mj/m 2 . 


/i 0 (nm) 

2 

4 

6 

8 

^rnB (/^n l) 

0.0102 

0.0236 

0.0547 

0.1382 
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For this case the linear theory predicts the inverse dependence of the dominant 
wave number with 7 and h 0 and hence with A (from equation (2.19) making 
P = 0). This shows the direct variation of dominant wavelength with these 
quantities. From nonlinear analysis it is found that bending moment can never 
overcome the destabilizing effect of LW force, so eventually rupture occurs in all 
cases. This implies that actual cell-support contact is possible for these cases. 
From fig. (2. 3) it is clear that the time of rupture (Tr) is minimum for a given A 
at some specific wavelength , which corresponds to I\ — K m - Tr increases with A 
(fig. 2. 3) for the dominant wave (K = K m ) as well as for the other wave numbers. 
This is due to the decrease of K m with increasing A (from linear theory) which 
again retards the growth rate of perturbation. Increase of A signifies an increase 
in the interfacial tension at constant mean film thickness (unlike h 0 , 7 has no 
contribution to the scaling parameter b for time), which makes the system more 
stable. The profiles at the time of rupture are shown in fig. (2.4) for different A at 
K m . The difference in the profiles is primarily due to the difference in K m which 
leads to different length scales of instability. The effect of different length scales 
of perturbation on the rupture profile for a given A , is shown in fig. (2.5) . 

2.6.2 CASE : 2 (P ^ 0). 

This is the case when both polar and apolar interactions are taken into account 
to evaluate the free energy. For all of the following cases, the cell-membrane is 
taken to be polar in nature (e.g. erethrocytes etc). According to the nature of 
the solid support this analysis is divided into two parts. For the first part' the 
substrate (the solid support) is apolar (e.g. teflon). So the polar interactions 
arise solely due to the polarity of the film-fluid and the cell-membrane. Typical 
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6 8 10 12 14 16 18 20 22 

A 

Fig. 2.3 Variation of time of rupture with wavelength in 
nondimensional form for different A. Apolar case. 





Fig. 2.4 Nondimensional profiles at the time of rupture for different 
A. Apolar case. The values of dominant wavelengths for the 
three curves are 8.26, 9.53 and 10.88 respectively. 





Fig. 2.5 Nondimensional profiles at the time of rupture for different 
wavelength. A=l. Apolar case. 





values of the components of the interfacial tension are taken such that S LW > 0 
('y BW > > 7 s W )i and S AB < 0 (77 > 7 c ~/4) (17). This leads to the system 

where AB interaction is destabilizing and LW interaction is stabilizing. For the 
next part the substrate is polar, thus at least one of the polar components of 
the interfacial tension is non-zero (e.g. agarose, cellulose acetate etc). In this case 
S LW < 0 (destabilizing effect) and S AB > 0 (stabilizing effect). As in the previous 
case the conditions for the different components of the interfacial tensions of the 
three members of the system are chosen according to the relations given in the 
appendix in (17). Here the first case resembles the region 4 and the second one, 
region 2 of the fig. 2 in (17). For the apolar support, the values of the spreading 
coefficients are: S LW = 2 mJ/m 2 , and S AB = —33 mJ/m 2 and those for the 
polar support are taken to be: S LW = —4 mJ / m 2 , and S AB = 25 mJ / m 2 . The 
properties of the fluid film is taken to be the same as those of water. 

The variation of the second derivative of the total free energy (<f>h) with respect 
to the film thickness is shown in fig. (2. 6). It is evident from the two curves that 
4>h is positive for the apolar case (fig.2.6a) beyond the film thickness 8.8 nm, 
whereas that for the polar case is below 8.1 nm (fig. 2. 6b). So these two curves 
are similar to the part in between hi and h 2 of the region 4 and region 2 of the 
fig.2 in (17), respectively. For both these cases a critical value of h — h c can be 
defined beyond which the competition between the two components of the free 
energy is more pronounced (17). The values of h c are estimated here at those 
points where the simulation results with the Fourier collocation technique show 
rupture (17) (though in the actual case there is indeed significant competition 
between the stabilizing and the destabilizing term(s), but no rupture is observed 
employing a efficient finite difference approach for discretization, which will be 
discussed in the next chapter). The results are obtained for the analysis in the 
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Following two cases with ho near or beyond the h c . as these cases are expected 
to be more interesting to study the additional effect of bending moment. 

The results from the linear stability study are shown in fig.(2.7). It deals 
with the variation of the critical wavelength A c at K = K c with the mean film 
thickness within the limit </>h < 0. For both these cases it is clear from the figure 
that the value of A c is much greater than the value of h 0 at any specific point 
on the curves, satisfying the criterion for long wave analysis which is given as 
A c /h 0 > 1. It is also obvious from this discussion that this criterion wall be 
automatically satisfied for the apolar case, since in that case the K c values will 
be smaller since P = 0 [2.18]. The general qualitative result from the nonlinear 
study is that in this case the actual cell-support contact never takes place. This 
finding is elaborated in the following subsections. 

• APOLAR SUBSTRATE: 

Fig. (2.8) shows the profiles at the time of rupture for different starting mean 
film thickness at constant interfacial tension. The apparent difference in the 
peak values of the profiles for different h 0 is again due to the difference in the 
length of the unit cell in which the instability evolves, owing to the decreasing 
K m for the different cases, with increasing h 0 . But the nondimensional time of 
equilibriation increases with increasing h 0 as shown in the figure (and hence the 
the dimensional time , t oc Th 0 \ [2.6]). This is because at higher ho the increase 
in the magnitude of P (oc h$) is much greater than that of A (oc hQ, thereby 
causing more destabilization. Fig. (2.9) shows the effect of 7 at constant ho- 
Here again the difference in the peak values can be physically resolved as earlier. 
The nondimensional time decreases with increasing 7 (also the dimensional time, 
t oc T [2.6]) as in the figure. This can be attributed to the stabilizing nature of 
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Fig. 2.7 Variation of critical wavelength with mean film thickness. 

7 = 0. The values of the spreading coefficients for 

the curves 1 and 2 are same as those in fig. 1.6. Polar case. 
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Fig. 2.8 Nondimensional equilibrium profiles of the film for different 
mean film thickness. Polar case with apolar support. 

7 = 0.001 mJ/m 2 , S LW = 2 mJ/m 2 and S AB = —33 mJ/rn 2 . 
The values of A, minimum nondimensional thickness at 
equilibrium and nondimensional dominant wavelengths for the 
four curves are given respectivaly as - Curve 1: 0.0244, 0.0199, 
4.36 ; Curve 2: 0.0144, 0.0246, 2.98 ; Curve 3: 0.007, 0.0346, 
2.27 ; Curve 4: 0.0023, 0.0587, 2.132. 
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Fig. 2.9 Nondimensional equilibrium profiles of the film for different 
values of interfacial tension. Polar case with apolar support. 

Mean film thickness = 5 nm. The spreading coefficients are 
same as in fig. 1.8. The values of A, minimum nondimensional 
thickness at equilibrium and nondimensional dominant 
wavelengths for the three curves are given as - Curve 1: 0.0, 0.019, 
2.986 ; Curve 2: 1.44, 0.020, 3.15 ; Curve 3: 14.44, 0.023, 4.78. 
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interfacial tension. However for both the above cases the dimensional time of 
equilibriation is in the order of micro- to milli- seconds which signifies the spon- 
taneous attainment of equilibrium. Again the nondimensional time of equilibri- 
ation for all of the above cases has the maximum order of 10 and the difference 
in the equilibiation time for the various cases is not significant and will not be 
observed in reality. So physically the stabilizing effect of the bending moment 
will not be that apparent. 

A specific study for the evolution of the cell-film interface with time for a 
typical case of 7 = 0.001 mJ/m 2 and h 0 = 5 nm at K — K m is shown in fig. (2. 10). 
It can be easily inferred from the figure that initially two drops are formed from 
the original perturbation which then coalesce to give one drop. Afterwards this 
drop starts thinning due to the destabilizing forces. Ultimately the stabilized 
equilibrium profile is attained (as evident from the fact that there is no change in 
it even after the lapse of long time) when the stabilizing components overcome the 
destabilizing ones. It is interesting to note that the time taken for the formation 
of the two drops from the original perturbation is much greater than that for 
equilibriation from the two drops. This is because at higher film thickness the 
body forces are insignificant in comparison to the interfacial tension which has 
slow effect on the process of evolution . But as soon as the small film thickness 
(at least at some points) is attained the body force terms become dominating, 
therefore dictate the ultimate fate of the film immediately. 

Among the six profiles shown in this figure, curves 1,2,4 and 6 are used to 
estimate the evolution of the complex nanostructures in three dimension, as given 
by the contour plots of fig. (2. 11). These structures are a combination of the holes, 
films, islands, and rims as anticipated in (19). Such structures are witnessed in 
reality in the stability of thin films for different systems (20,21). 
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Fig. 2.10 Evolution of film profile with time. Mean film thickness = 5 nm. 
Interfacial tension = 0.001 mJ/m 2 . Spreading coefficients are 
same as in fig. 1.8. A — 0.0144. The nondimensional dominant 
wavelength = 2.98. Minimum nondimensional film thickness at 
equilibrium = 0.0246. Polar case with apolar support. 
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than 0.2 



H = 0.2 -0.9 


H = 0.9- 1.0 


H = 1.0- 1.1 


H = 1.1 -1.5 


H greater 
than 1.5 


Fig. 2.11 The contour plots for the 3D configurations of the selected 

curves of fig. 1.10 drawn for 4A m : 4A m . The fig (A), (B), (C), 
and (D) respectively correspond to the curves 1,2,4, and 6 of 
fig. 1.10. 







They are supposed to be important in multilayer adsorption and heterogeneous 
nucleation (19). Since the consideration of nonlinear dynamics of the films for 
two dimensional perturbation in a large domain is difficult with the present com- 
putational facilites, the method of superposition of waves with random phase 
shifts and orientation (22) is adopted to generate these structures. At any point 
in xy plane, the net contribution to the film thickness (Ht) by M number of 
waves of constant wavelength, which are randomly oriented and phase shifted, 
can be written as: 


M 

H t = '%2f[K (x cos 0j + y sin 6j) + <j>j] (2.38) 

j=i 

where 9j is the orientation of the jth wave, <f>j is its phase shift and K is the 
■wave number which is constant for all of these waves. The contribution of each 
of these waves can be obtained from the information available for the x vs H 
(simulated film thickness) profile using its periodic property and suitable inter- 
polation technique. Then the film thickness at the said point can be estimated 
by dividing Hj with a suitable scaling factor (SF) . SF can be calculated from 
the mass conservation principle as given by: 

SF = S§- (2.39) 

Here Ht is the mean value of Ht, calculated by using Ht values for all the 
points considered in the xy plane. H is the mean value of the film thickness JT, 
obtained from the simulation. 

The four snap shots at different time show how the small holes are devel- 
oped and then propagated to produce larger ones. To describe it more clearly, 
fig.(2.1lA) shows the condition at initial perturbation. Then the second one 
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(fig.2.1 IB) shows the intermediate evolution containing larger variation in the 
local film thickness, as is evident from the increase in the number of shades than 
in the previous one. The third figure (fig. 2 . 1 1 C) corresponds to the two drops 
configuration in two dimension. It shows larger holes (increased depth and width) 
because of the appearance of two more shades and the reduction of the area oc- 
cupied by the initial two shades. Fig. ( 2.1 ID), corresponding to the equilibrium 
profile, gives these holes with much larger dimensions since the effect on the 
shades, as discussed for the third figure, is more prominent here. 

The last part of the results for this case, as shown in fig. (2. 12), deals with 
the evolution of a random perturbation for the same values of h 0 and 7 as in the 
previous result. The initial and final profiles are shown for a length proportional 
to that of three unit cells. The stabilizing nature of the system is more confirmed 
in this case. 

• POLAR SUBSTRATE: 

Only one parameter study is shown for this subgroup as the other ones are 
almost same as those in the previous case. Fig. (2. 13) shows the variation of 
the equilibrium profiles for different ho at constant interfacial tension. The only 
interesting thing can be discussed here is that the thin film structure in the 
ultimate profiles are almost in the order of the original mean thickness. The 
other features can be explained similarly as in the previous case. 

2.7 CONCLUSION AND RECOMMENDATIONS 

It can be concluded that, for low interfacial tension cases bending moment has an 
important role in the film stability and hence has to be incorporated to study the 
cell-support adhesion phenomena. With bending moment and interfacial tension 
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Fig. 2.12 The initial and final profiles for a random purturbation. The all 
parameter values are same as in fig. 1.10. The figure is drawn 
for 3A m . Minimum nondimensional film thickness = 0.025. 
Nondimensional time of equilibriation = 10. Polar case with 
apolar support. 
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11 

Fig. 2.13 Nondimensional equilibrium profiles for different mean film 

thickness. Polar case with polar support. Interfacial tension is 
0.001 mJ/m 2 . S LW = —AmJ/m 2 and S AB = 25mJ/m 2 . The values 
of A, nondimensional minimum film thickness and nondimensional 
dominant wavelengths are given as - Curve 1: 4.08, 0.061, 18.34 ; 
Curve 2: 3.31, 0.072, 16.67 ; Curve 3: 2.61, 0.083, 15.06 ; Curve 4: 
2.0, 0.098, 13.54. 
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together, the qualitative results do not change much from the case of instability 
study for the gas-liquid film interface. But quantitatively for different values of 
interfacial tension, the results have shown a large varitation, specially for the 
length scale of instability and the ultimate structure of the cell-film interface at 
the time of equilibriation or rupture. This is expected since interfacial tension 
has stabilizing effect on the film. The time scale variations can be neglected as 
they are micro seconds phenomena. So these cell membranes evolve in a larger 
length scale to form a more stabilized ultimate configuration. 

From the observations of this chapter it can be proposed that, to study the 
physics of this type of systems, one may include other effects like Marangoni 
effect, the effect of presence of some specific molecules on the membrane surface 
etc., together with the effect of bending moment to get a more realistic result, 
according to the nature of the study. One can also include the effect of Born re- 
pulsion (23), but it may not be useful as the system is originally stable for most 
of the cases discussed above. Another important thing to be noted is that for low 
surface tension and highly curved surface structural effect, like the existence of 
a structural double layer can play an important role in the interfacial instability. 
Further work may be carried out in this respect with the complete stability anal- 
ysis of highly curved and charged films (5). Since the resulting profiles developed 
in case of the bending moments are quite interesting and varying for different 
cases, in order to visualise their three dimensional form actual three dimensional 
analysis may be done, instead of the statistical estimation of those structures 
presented here. 
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Chapter 3 


EFFECT OF RECEPTOR 
AGGREGATION ON THE 
MEMBRANE SURFACE IN 
THE PROCESS OF 
CELL-SUBSTRATE 
ADHESION 

3.1 INTRODUCTION 

In the study of cell adhesion on a solid support another important aspect is the 
presence of specific molecules on the surface of the cell-membrane called recep- 
tors. It has been proposed by several workers (2,9,24) that at the moleculer 
level cell adhesion to the support may be mediated by a wide variety of adhesion 
molecules, termed as receptors, interacting with a broad spectrum of exogenous 
surfaces. They are also thought to be possessing adhesive properties with respect 
to some specific sites on the solid support. In the present work this particular 
concept is viewed in the light of modified specific interactions between the mem- 
brane and the support due to the presence of receptors. The receptors try to 
aggregate at some points or other on the membrane under the influence of both 
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ordinary diffusion and convection caused by local membrane curvature. This 
aggregation in some cases may change the original interaction properties of the 
bare membrane, resulting in a different type of system behaviour. Here a model is 
proposed to modify the interaction potential and thus the corresponding terms in 
the governing equations of the system (evolution equation for the membrane-film 
interface and mass balance equation for the receptors). Using proper techniques, 
numerical experiments are carried out for some specific cases to validate the above 
proposition. 

3.2 MODELING OF INTERACTION POTEN- 
TIAL IN PRESENCE OF RECEPTORS 

The basic tenet for this modeling is that the receptors can aggregate and change 

the potential of the membrane locally. At some points the receptor concentration 

(r) becomes maximum. The other limit is the value for the bare membrane at T = 

0. The behaviour can be thought of being linear initially, and after some certain 

value of T, becomes asymptotic when the modified potential of the membrane 

> 

is close to that of the receptors. Mathematically this can be expressed by the 
following equation: 


AG = A G„ + (Agr c + A r Gm) - (3.1) 

where A G m is the free energy change for the bare membrane and A G r , that 
for the receptors. C is the modeling parameter which dictates the value of Toe. The 
variable AG is the modified free energy change for the system which is equivalent 
to the spreading coefficient described in Chapter 1. Similarly the values of A G m 
and A G r can be evaluated from the corresponding data of interfacial tensions, 
as in Chapter 1. The above equation [3.1] can be used for both the LW and AB 


46 



components of the free energy. 


3.3 FORMULATION 

The basic system discussed here is same as that in Chapter 1, except that the cell 
membrane now contains some receptor molecules which can interact with some 
binding sites on the solid support (fig. 3.1) . The extra body force coming into 
play here (other than those stated in the previous chapter) is the attractive force 
between the receptor molecules and the binding sites. This force modifies the 
general definition of the body force term in the Navier-Stokes equation. 

Starting from a reference state with an undisturbed film of thickness h 0 and 
receptor concentration of To (number of receptor molecules per unit area), the 
dynamics of the system can be modelled by using N-S equation and modified 
continuity equation to predict the variation of F. A long wave analysis is done 
as before to get the evolution equation as follows: 


ht + 


1 

12/1 



(3.2) 


where the meaning of each term is the same as the Chapter 1. However, unlike 
the (j) in Chapter 1 , here <j> is a function of both h and T . So (f) can be expressed 
as follows: 


, 2 A G LW <P 0 A G AB (d 0 — h,' 

<t>= ^ 5 ex P 


h 3 


lo 


k 


(3.3) 
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BINDING SITES 


Fig. 3.1 Configuration of the system 


48 



where A G LW and A G AB are functions of F, given, in general, by equation 
[3.1]. Considering this functional relation of <f>, the above evolution equation [3.2] 
can be modified as: 


ht ~t~ 22^ hxxx h ( (j>h hg 4 " ^rF^)!^ — 

The mass balance (convection-diffusion) equation for the 
follows 


0 (3.4) 

receptor is given as 


T t = D S T XX + (Gr) x J (3.5) 

where Ds is the mass diffusion coefficient of the receptor molecules, a; is 
the area occupied by each receptor molecule (surface coverage), / is the friction 
factor, and G> is the total free energy for the receptors. / is defined as the ratio 
of KT/Ds where K is the Boltzman constant and T is absolute temperature of 
the system. G T is given by the following equation: 


G r 


A G^ w d 2 0 
fi 2 


+ AG^ s exp 



(3.6) 


where A G r can be defined as the change in free energy of the receptors. [G r ) x 
can be substituted by ( G r ) h h x . 
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The first term in equation [3.5] is the diffusive term and the second one is the 
convective term. The driving force per molecule is given by — ( G T ) X . So the total 
force for all the molecules can be obtained as —Fa,- ( G r ) x . The total force is then 
divided by the friction factor to obtain the velocity. In this way the convective 
nature of the second term can be explained. 

• NON-DIMENSIONALISATION OF THE EQUATIONS 

The above equations are non-dimensionalised to give 


H T - (A OS') {* (l + ^Py) - *** exp (?) (> + ^j) } * 

+ [spn (ACT) { - 0 ? exp ( V) } Fx + iTfe] = 0 


Tt = ^6 


A t T xx - sgnAG, w { T 


W-L- 

{H 3 


^• 5 exp (-J-) | H x 


The following non-dimensional parameters are defined: 




(3.8) 


// — SSL 1—lSL 
a ~ h 0 5 L ~ h 0 

where v = ~ . The definitions of a and b are as follows: 

p 


(3.9) 


a = 


6cP 0 |AGJf| L (6<P„|AGr|)' 

^ — ’ 6= iwr 


Other dimensionless constants are given by: 


(3.10) 


p c _ LW _ (AG, W A r ,AB _ ( /".G; ,B \ . 

C ~T 0 ’ a -[a ^- 1 ]’ 0 ~ \AG^ B V (3 ' U) 
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A _ hj £G£ A CG 


LW 


3r 0 


r 0 hl(AG? B -AG£B) A _ fDs hl 
~~ 6Cl 0 dlAGW ’ A 4 - 2a,dg|AG^|’ 


(3.12) 


. _ , _ 4ain\AG BW 

yls ~ 2J 0 i2AGi^’ 716 ~ /io/|AG^| 

The six constants Ai to A& [3.12] help to understand the behaviour of the 
equations [3.7] and [3.8]. A\ shows the relative effect of the LW and AB forces 
for the bare membrane. Ai and A 3 govern the effect of the receptors, A 3 describes 
the polar contribution of the receptors. A 4 shows the relative effect of the diffusive 
part compared to the convective part in the mass balance and A 5 has similar effect 
as that of A\ in the convective part. A& has the important role of predicting the 
time scale of aggregation of the receptor molecules. The other two constants 
G LW and G AB couple the potential field of the receptors with that of the bare 
membrane. In all, the term in the first braces in equation [3.7] has stabilizing 
or destabilizing effect depending on the values of sgn(AG^ y ) , A\, G m and 
G AB . The first term in the second braces shows the effect of stabilization or 
destabilization depending on the values of sgn^AG^j, A 2 and A 3 . These two 
terms are coupled by the parameters H and T and their respective derivatives. 
In the mass balance equation [3.8] the first term is diffusive and the one in the 
braces is convective. 


3.4 LINEAR STABILITY ANALYSIS 

As in Chapter 1, the two governing equations for the overall system [3.7] and 
[3.8] can be linearised around the mean values of h and F (i.e. h 0 and To). In 
that case the solution of the nondimensional evolution equatioq can be written 
as follows: 

CENTRAL library 

». I.T., KANPUR 

zz* A.nmi 
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H — 1 + Eh exp (iK X + fiT) (3.13) 

and that for the nondimensional mass balance equation is of the form: 


r — 1 H“ Er exp (iK X -f- fiX*) 


(3.14) 


where Eh and Er are the nondimensional amplitudes of the perturbations 
applied to the system given by: 


= ^ (3.15) 

flo to 

where th and ep are the dimensional amplitude of perturbations. K is the 
wave number and Cl is the growth rate of the perturbations in nondimensional 
form. 

Now using these solutions, equations [3.7] and [3.8] can be linearised to give 
the expression for Er and Cl in terms of the nondimensional constants (shown 
before) and the wave number K as follows: 

A 6 sgn (A G™) (l - exp (^)) E h K 2 
Er ~ fi + A 4 A 6 I<i 

n 2 + Cl n + c 2 = o (3.17) 

where Ci and C 2 are given by: 



ci = K 2 


K 2 + sgn (A G L m W ) 1 + 


g lw 

C + 1 


A 1 exp 



+ A 4 Aq 


(3.18) 
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c 2 = A,A e K< [IP + sgn ( AG™) { (l + §^) - A, exp (^) (l + §^) }] 

+ S S™ (A G™) sgn (AG™) (l - 4. exp (AJjj (4 2 - 4 3 exp (- 4 1 )) 

(3.19) 

The expression for the nondimensional amplitude of receptor perturbation 
implies that for the cases with attractive van der Walls force, the evolution of 
receptor concentration will be just the opposite of that for the film thickness. 
So larger concentration will be found at the smaller film thickness, which will 
significantly affect the free energy change at those points. The has two roots 
of which only the positive one is responsible for the instability. The value of K 
corresponding to the maximum value of gives an estimate for the wave number 
of the fastest growing wave. The set of values of 0 and K for this case can be 
used to obtain the linear time of rupture as shown in Chapter 1 . 

3.5 NON-LINEAR ANALYSIS 

The linear theory only gives the estimate of the initial (upto which the system is 
fairly linear) behaviour for any system. So the better way to study the physics 
of the system is by using the nonlinear approach. The observations are broadly 
divided into two parts. In one case the system containing bare membrane along 
with the film and the support is originally apolar, and in the other case it is 
polar. 

NUMERICAL METHOD 
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Here the method of solution of the two coupled partial differential equations is 
almost same as that in Chapter 1, except that for the purpose of discretization the 
approach of finite difference is adopted instead of Fourier collocation technique. 
The reason of this change is that the system becomes highly nonlinear after 
some time at some specific locations. This needs high grid density for these 
locations. Use of Fourier collocation for such high number of grid points is very 
time consuming and at the same time gives erroneous result. 

The concept of finite difference used here is given below: 

Yi + i - Y { _ i 

PH = J AX ’ (3-20) 

where Y is any fuction value and Yx is derivative of Y with respect to X at 
the 2 th grid point. The necessary function values at the i + | and i — \ grid points 
are interpolated using suitable technique. Like in Chapter 1, here also the inital 
condition is taken from the solution of the linear stability analysis at T = 0 and 
periodic conditions are used as the boundary conditions. 

3.6 RESULTS AND DISCUSSION 

To obtain the results for different cases considered so far, the specific parameter 
values for the system are taken to be of the same order as those used in Chapter 1. 
The free energy values for the receptors for both LW and AB are also considered 
to be of the same order as those of bare membranes. Here two specific cases are 
considered regarding the bare membrane to make a comparative study : firstly, 
the apolar case with strong apolar attraction of van der Waals force and no 
polar part, and secondly, the more interesting one, the polar case where the van 
der Waals attraction together with the repulsive polar components are present. 
The actual values of the free energy change (spreading coefficients) in these two 
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cases for membrane, used in computaion, are shown in table (3.1). The receptor 
potential considered for these two cases are same, given by A G^ w = —4 mJ /m 2 
and A G^ B = —45 mJ/m 2 , that is, both the components are destabilizing in 
nature. The values of the parameters (7) other than those used in Chapter 1 are 
summarised in table (3.2). The values of a,- and C are varied to get different set 
of results for each of the above two cases. 

At first linear analysis is done for both the cases to get the values of K m for 
them. The variation of Cl with K and ho for the two cases is shown in fig. (3.2). 
For the apolar case, the disturbance is maximum at small values of ho. On the 
other hand, for the polar case this is observed for h 0 around 9 nm . To draw a 
better comparision of the results of polar and apolar cases, h 0 is assigned a value 
of 10 nm (at which the instability for the polar case is the least among the values 
shown in fig. (3.2)) in both the cases. 

For the nonlinear results, the grid convergence study is done at first. The 
comparative results for the polar case with 400 and 800 grid points in the full 
wavelength are shown in fig.(3.3). It shows the ultimate h and T profiles. The 
results for the two sets of grid points are fairly similar except for the values of 
r max • The time taken to reach this condition is also almost same for these two 
cases. Though the profiles for the 800 grid points are smoother, but this requires 
approximately triple the computational time. So for each of the two cases stated 
above, only one set of results will be presented with 800 grid points later. 

The results of the nonlinear simulation for the apolar case is shown in figs. (3.4) 
and (3.5) for one set of values of a; and C. The LW component of free energy for 
the bare membrane is taken here to be double of that for the receptors, to make 
the system highly apolar in nature. Fig. (3.4) shows the variation of the minimum 
film thickness and maximum receptor concentration with time and fig. (3.5) shows 
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Table 3.1: The values of different components of free energy in mJ/m 2 for the 
bare membrane 



LW Component 

AB Component 

Apolar Case 

-8 

0 

Polar Case 

-4 

25 


Table 3.2: Estimated parameter values 


Parameter 

Best estimates 

Absolute Temperature 

298 K 

Viscocity of film-fluid 

0.8 cp 

Surface diffusion coefficient 

10 -13 m 2 /s 

Surfactant concentration 

10 12 — 10 15 molecules/m 2 

Surface Coverage of surfactant 

10-1000 ( nm ) 2 
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(a) (b) 

Fig. 3. 2 Variation of dimensional growth rate with dimensional 

wave number for different values of h 0 ■ For receptor 
A G LW = —4 mJ/m 2 and A G AB = —45 mJ/m 2 . 

(a) Apolar case, i.e., for bare membrane /S.G LW = — SrnJ/m 2 
and A G AB = 0. (b) Polar case, i.e., for bare membrane 
A G LW = — \mJjm 2 and A G AB = 25 mJ/m 2 . 

Values of a,- and C are 10 nm 2 and 2, respectively. For 
curves 1 to 5, ho values are 8, 8.5, 9, 9.5, 10 nm, respectively. 
Interfacial tension is the stabilising factor. 
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0 


Fig. 3.3 Nondimensional profiles at the time of rupture for (a) 400 
grid points, and (b) 800 grid points. Polar case with h 0 
= 10 nm. Other parameter values are same as in polar case 
of fig.3.2.. Interfacial tension is the stabilising factor. 
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T 


Fig. 3.4 Variation of nondimensional minimum film thickness and 

maximum receptor concentration with nondimensional time. 
Apolar case with h 0 = 10 nm. Other parameter values 
are same as in apolar case of fig.3.2 .Interfacial tension is • 
the stabilising factor. 
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0 1 2 3 4 5 6 

(a) 


3 3.05 3.1 3.15 3.2 3.25 3.3 

*n 

(b) 


Fig. 3.5 Nondimensional profiles at and just before rupture. 

(a) just before rupture and (b) at the time of rupture. 
Apolar case with ho = 10 nm. Other parameter 
values are same as in apolar case of fig. 3. 2. Interfacial 
tension is the stabilising factor. 
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the profile of the film thickness and receptor-concentration just before and at the 
time of rupture. It is clear from these figures that the receptor aggregation has no 
effect on the film stability. For both the cases with and without receptors, the film 
ruptures. However, quantitavely, just before the rupture, due to high gradient 
of the film thickness with the length, near the point of rupture the receptor 
concentration increases (fig (3.5a)). This accelerates the process of rupture and 
the receptor concentration shoots up to a high value compared to the mean 
concentration at and around the point of rupture. So on the whole, in this case 
the receptors cannot affect the bare membrane potential so much as to change 
the process of cell adhesion (both the dynamics and equilibrium solution) on the 
support. The experiments for the other combinations of a; and C show that they 
can only change the maximum receptor concentration at the time of rupture and 
the time of rupture itself. 

The results of the nonlinear studies for the polar case is interesting. Here the 
LW component of the free energy for both the bare membrane and the receptors 
is taken to be the same in magnitude and sign to decouple the effect of receptors 
on the AB component of the free energy for the system. As shown by some 
workers, the film always attains an equilibrium profile, made up with ultrathin 
film and drop(s), in those cases where LW and AB components of free energy are 
opposite in nature (17,19). The same observation may be extended for the case 
of bare membrane system with same type of potential, as is found in the previous 
chapter. But the results with receptors show that the film ruptures and hence the 
actual contact of cell and support becomes possible, which is qualitatively just the 
opposite. This can be rationalised from the view point of receptor aggregation at 
some points of the membrane during its evolution. A nondimensional free energy 
term can be defined here by dividing equation [3.1] by A G m . This can be applied 
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for both the LW and AB components of free energy. Since here only the effect 
on AB component is taken into account, it is represented in nondimensional form 
by: 


A G AB = 1 + £ — (3.21) 

§ + l 

It can be easily identified with the second expression in the braces of the first 
term of equation [3.7]. The study of this term during the process of evolution will 
help in understanding the force field acting on the system at different instant. 

With the help of fig. (3.6), (results with same set of values of a,- and C as 

that in the apolar case, shown in the nondimensional forms), the process of 

evolution for a typical case can be explained. Initially when there is no significant 

receptor aggregation, the membrane behaviour is similar to that of the bare 

membrane. The steady profile is attained after some time of the initial perturbed 

condition which is almost unchanged for a long time, during which the receptor 

aggregation takes place. This is shown by the flat region of the H x curve in 

fig.(3.6d). The curves 2 show the condition at an intermediate time in between 

initial perturbation (curves 1) and the steady state profile (curves 3). The curves 

3 are drawn at a specific instant during this period of time. .Upto this point 

there is a slow increase in the maximum value of T as shown by the r max curve 

of fig.(3.6d). Then T attains a high value at some point(s) on the membrane 

surface which further facilitates the evolution, as shown by the curves 4. The 

ab 

interesting difference between the AG profiles of curves 3 and 4 is that, for the 
latter the free energy attains negative values at some points on the membrane, 
which signals the recommencement of the instability. Now both h and T affect 
each other. The minimum film thickness decreases, while the maximum receptor 
l concentration increases at that point, as shown in the curves 5 of fig.(3.6c) and 
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(c) (d) 


Fig. 3.6 Evolution of nondimensional profiles (a) to (c) with nondimensional 
time, (d) Variation of the nondimensional minimum film thickness 
and maximum receptor concentration with nondimensional time. 
Polar case with ho = 10 nm. Other parameter values are same 
as in polar case of fig.3.2. The nondimensional times for the curves 
1 to 6 are T = 0, 8, 225, 241.4, 241.8, 241.8009, respectively. 
Interfacial tension is the stabilising factor. 
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by the rightmost parts of both of the curves in fig.(3.6d). Now the potential of 
the system is considerably affected by the receptor concentration which makes 
the change in AG (here only A G AB ) values at some points highly negative (see 
curve 5). This ultimately leads the film to rupture, together with high receptor 
concentration at the rupture point (the curves 6). The curves obtained from the 
simulation result are not smooth due to less number of data points. A smoothed 
version of the curves 6 of fig. (3.6c) is shown in fig. (3.7). 

As in the previous chapter a 3D structure is generated for the above case using 
the principle of random wave superposition. Three specific instants are chosen, 
namely, at the time when quasi-equilibrium state is attained when A G AB is non 
zero for all points, then at an instant when the finger like structure is already 
developed and A G AB is negative at those points, and at an instant just before 
rupture. Six snap shots are presented in fig. (3.8), the left column corresponding 
to the film thickness and the right one to the receptor concentration at the same 
instant. In fig.(3.8B) there is very small variation of receptor concentration which 
shows the quasi equilibrium state. In fig.(3.8D) the receptor concentration varies 
largely which shows that the process of aggregation starts. Also in fig.(3.8C) 
(corresponding to fig.(3.8D)) the film thickness attains lower values and higher 
values at some points or other which shows that the process of evolution again 
starts. Fig.(3.8F) shows black points spreaded in throughout region in a random 
fashion which means maximum receptor concentration at these points and hence 
rupture will occur at these locations. This is also supported by fig.(3.8E) which 
shows white regions (minimum thickness) corresponding to those black points. 

A wide range of results is obtained for this case when different combinations of 
a t - and C are tried. As seen from equation [3.1], as C increases, AG (here A G AB ) 
becomes more insensitive to the variation of F. So, AG needs more receptor 
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T] 

Fig. 3.7 Curves 6 of fig. 3. 6 after smoothing. 
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Contour plots of the 3D profiles, estimated by using superposition 
of randomly oriented and phase shifted waves, for the case same 
as in fig.3.6. The left column for H and right one for T. 

(A) and (B) are at the quasi steady state. (C) and (D) are at an 
instant when receptor aggregation started. (E) and (F) are just 
before rupture. The key for the grey scale is given in fig. 3. 9. 













concentration to change from positive to negative value. Increase in C causes 
more aggregation and also prolongs the time of rupture, as is clearly seen from 
table (3.3). But there is a limit to which C can be increased (Coe) for a given 
value of ai, because after that receptors can no longer change A G significantly 
( C >> T), and so another equilibrium solution is reached. This can be avoided if 
the amplitude of the perturbation for the receptor concentration (£r) is increased 
(by increasing the value of a,) so that considerable change of AG can occur even 
at high C . For example, for a,- = 160 (nm) 2 at C — 2.7 there is no rupture, 
whereas for the same C value at a,- = 250 (nm) 2 there is rupture at T = 68 
with r ma x = 5.5. Again, for a given C, as a,- increases the receptor aggregation 
as well as the time required for rupture increases (table (3.3)). This is due to 
the fact that a,- increases the value of A 6 [3.12], which dictates the time scale of 
receptor- aggregation [3.8]. For a t also, there exists a saturation value after which 
no significant change in rupture time or the receptor aggregation is possible. 

Some interesting results are shown in fig. (3. 10) to demonstrate the combined 
effect of a{ and C on the receptor aggregation and film instability. Their values 
are increased for the curves 1 to 3 shown in this figure. r max is found to increase 
gradually in all the cases, but the time of rupture (and hence the flat region ) 
decreases from curve 1 to 2 and then increases again from curve 2 to 3. The 
decrease in flat region in curve 2 shows that for that value of a,- higher value of 
C can be used before C K is reached. The large increase in both rupture time and 
maximum receptor concentration at that instant in curve 3 compared to that in 
curve 1 indicates the mutual reinforcing effect of these two variables. The region 
of the profiles where considerable changes take place near and at the time of 
rupture is shown for these two curves in fig. (3. 11). The larger negative value of 
AG' for curve 3 compared to that for curve 1 causes more receptor aggregation. 


68 



Table 3.3: Results from simulations with different values of a,- and C 


a.i (nm) 2 




C 

r 

1 max 

Tr 

160 

2 

3.38 

5.3 

160 

2.2 

3.67 

10.28 

160 

2.6 

4.28 

893.56 

170 

2.6 

4.38 

482.65 

200 

2.6 

4.72 

68.473 

250 

2.6 

5.3 

32.7 





Fig. 3.10 Variation of nondimensional minimum film thickness 

and maximum receptor concentration with nondimensional 
time for different parameter values. Polar case with h 0 
= 10 nm. The LW and AB components of A G for bare 
membrane and receptor are same as in fig.3.6.. Values of a,- 
and C for the curves 1 to 3 are respectively - Curve 1: 10 nm 2 
and 2; Curve 2: 250 nm 2 and 2.7; Curve 3: 500 nm 2 and 3.3. 
Interfacial tension is the stabilising factor. 
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• BENDING RESISTANCE AS STABILISING FORCE 


So far, interfacial tension is considered to provide the stabilising effect to 
the system. As discussed in the previous chapter, this is valid for large value 
of interfacial tension at the membrane-film interface. But for small value 
of this parameter, the effect of bending resistance during the deformation 
of the cell has also to be taken into account. The behaviour of the system 
is now studied by considering the bending resistance as the only stabilising 
factor. The eq. [3.2] can be modified in this case as follows: 


h t + 77— (—Bh 3 h xxx - h 3 <f> x ) = 0 
12 n v 

The nondimensional form of the eq. [3.22] can be written as: 


(3.22) 


H t 


+ 


H x 


( AG “) {* (* + ypy) - ***»«• M l 1 + pj 

sgn (AGJf ) { exp (^)}r* - iPH xxxxx ^ = 0 

X (3.23) 


where the value of space scaling parameter a will be given as: 


(3.24) 


The other parameters and equations remain the same. From linear stability 
analysis the dispersion equation for this case can be wriiten similar to that of 
[3.17]. But the expressions of C\ and C2 will be as follows: 


C] = A' 2 


K 1 + sgn (A Gi w ) 


1 + 


G 


LW 


C + 1 


Aiexf (izl)U + G 


;AB 


l 


C + l 


+ A4A6 


(3.25) 
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<* = A 4 A 6 K 4 + s g n (AGST) { (l + $£) - A exp (fcl) (l + £) }] 

+ (l«7^ n ( AG ™ W ) s S n ( AG ^) { l - 4> exp (^i)) (x4 2 - Az exp (*=*)) 

(3.26) 

With the help of the above modified equations, simulation is carried out for 
perturbation with dominant wavelength. The spreading coefficients are taken to 
be same as in polar case with interfacial tension. The values of a t - and C are taken 
as 10 (nm) 2 and 2 respectively. 400 grid points are used, all other variables being 
same as in the previous case with interfacial tension. The results are shown in 
fig. (3. 12). In this case also the receptor aggregation causes the rupture of the film 
and hence cell-substrate contact is made possible. Fig. (3. 12a) shows the variation 
of H m in and F max with T which is similar to the curve (d) in fig.(3.6). But from 
fig. (3. 12b) it is clear that in this case more aggregation of the receptor is needed 
to cause the rupture compared to the case with interfacial tension. This is due to 
greater stabilising effect provided by the bending resistance than the interfacial 
tension. 

3.7 CONCLUSION AND RECOMMENDATIONS 

Receptors are seen to have a significant effect in case of cell-support adhesion as 
is clear from the results presented above. For some specific cases their effect may 
be negligible (as for the apolar case), but for some, they may change the 
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Fig. 3.12 (a) Variation of nondimensional minimum film thickness and 

maximum receptor concentration with nondimensional time, 
(b) Nondimensional profiles at the time of rupture around the 
position where rupture takes place. Both are for polar case, 
all parameter values being same as in fig. 3.6. and bending 
is the stabilising factor. 
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entire qualitative behaviour of the system (the polar case). The concentration 
and surface coverage of the receptors are important in these cases. The model 
adopted here is found to be adequate to represent the actual physical phenomena 
and flexible enough to account for different cases with largely varying system 
parameters. Further works can be done in this respect using different approach 
to describe the modified interactions. Several other parameter cases can be tested, 
even with this model. The effect of bending moment, Marangoni effect, variation 
of the surface diffusion coefficient etc., can be included in these studies. Moreover, 
this approach can be generalised for the cell-cell interactions. 
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